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The flux penetration near a semicircular indentation at the edge of a thin superconducting strip 
placed in a transverse magnetic field is investigated. The flux front distortion due to the indentation 
is calculated numerically by solving the Maxwell equations with a highly nonlinear E(j) law. We 
find that the excess penetration, A, can be significantly (~ 50 %) larger than the indentation radius 
ro, in contrast to a bulk supercondutor in the critical state where A = ro. It is also shown that the 
flux creep tends to smoothen the flux front, i.e. reduce A. The results are in very good agreement 
with magneto-optical studies of flux penetration into an YBa2Cu30x film having an edge defect. 



I. INTRODUCTION 

Magnetic field penetrates type-II superconductors as 
a set of quantized flux lines - vortices. Macroscopically, 
the vortex matter can be considered as a "flux liquid". 
An important feature of this matter is pinning of vortices 
leading to zero electrical resistance at zero temperature. 
The pinning results in a non-uniform distribution of mag- 
netic flux forming a critical state. The critical state de- 
termines the macroscopic properties, e.g. the maximum 
current density and magnetic susceptibility, that are im- 
portant for applications. According to the critical state 
model, 1 at any point of the sample the local value of the 
electrical current density is equal to its critical value, j c , 
for a given magnetic field and temperature. 

An interesting property of the critical state is that lo- 
cal material defects affect the field and current distri- 
butions on a global scale. For example, even a small 
non-superconducting cavity or an edge indentation cre- 
ate sample-spanning discontinuity lines where the current 
flow direction changes abruptly. 2 At a non-zero tempera- 
ture, the critical state is relaxed due to flux creep that is 
conventionally described by a highly nonlinear current- 
voltage curve, E oc j n where n ^> 1 and E is electric 
field. Nevertheless, the same tendency persists: a small 
cavity of size i in a bulk superconductor perturbs the 
field distributions on a much larger scale of ~ n£. 3 Many 
applications of superconductors are based on thin films 
where this tendency must be even stronger since the rela- 
tion between the magnetic field and current is nonlocal^ 
Usually this leads to poorer performance of supercon- 
ducting devices whose global properties are deteriorated 
by numerous natural defects blocking the current flow. 
However, the same tendency can help control the flux 
motion on a global scale by patterning the superconduc- 
tor with arrays of small holes designed, e.g., to guide the 
flux in a particular direction^ 

Surprisingly, a quantitative understanding on how a 
single local defect affects the flux penetration into a su- 
perconducting film is still rather poor. Even the sim- 
ple case of an infinitely long thin strip with a semicir- 
cular edge-indentation is not solved. For a bulk super- 



conductor in the critical state such an indentation cre- 
ates an excess flux penetration exactly equal to the in- 
dentation radius #2 However the nonlocal electrodynam- 
ics in thin films and hence the presence of Meissner cur- 
rents in the flux-free regions make the picture much more 
complicated^ It has been observed using magneto-optic 
imaging that the excess flux penetration in films can sig- 
nificantly exceed the size of the indentation it originates 
from^ The physical mechanism behind this enhancement 
is however not yet understood. It could be related to the 
effect of thin-film geometry, to the flux creep or to ther- 
mal instabilities nucleated at the indentation. 

This work aims to clarify this question by presenting 
a detailed study of flux penetration into a strip with a 
semicircular edge-indentation in the flux creep regime. 
We determine how the excess penetration depends on the 
size of the indentation, the applied magnetic field and the 
creep exponent n. 



II. MODEL 

Consider a thin superconducting strip of thickness d 
placed in a transverse magnetic field. The strip is in- 
finite in the ^-direction, has the width 2w ^> d in the 
x-direction and a semicircular indentation with radius ro 
at the edge, see Fig. [TJ The flux dynamics in the creep 
regime is conventionally described using a local relation 
between electric field E and current density j^ i 1Q i n 



E = pj, 

with a highly nonlinear resistivity 
P = Po (j/jcT~ 



(1) 



(2) 



which does not explicitly depend on the magnetic induc- 
tion B. Here po is a constant, j c is the critical current 
density, while n is the creep exponent, n ^> 1. This expo- 
nent can be related to the activation energy U for ther- 
mal depinning as n ~ U jkT . Hence, large n means small 
creep, and the Bean critical state model 1 is regained in 
the limit n — > oo. 
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FIG. 1: A superconducting strip of width 2w with a semicir- 
cular indentation of radius ro in transverse field H a . 



For numerical simulations of flux penetration 
into the strip we use the formalism developed by 
Brandt^ii^^iiJ^il that can be applied to thin 
type-II superconductors of various shapes. For a thin 
superconductor, it is appropriate to look at length 
scales larger than the thickness <i, and introduce a sheet 
current 



rd/2 

J(r) = / dzj(r,z), 

J-d/2 



where r = (x, y) are in-plane coordinates. Due to the 
current conservation, V • J = 0, the sheet current can be 
expressed through a scalar function g(r) as 



J = V x zg 



(3) 



where g has the interpretation of the local 
magnetization. 9 Substituting the current from Eq. ([3]) 
into the Biot-Savart law one arrives at a non-local 
relation between B z and g, 

B z (r, z) = ^H a + f dV Q(r, r', z) g(r') . (4) 

J A 

Here H a is the applied field and A is the sample area. 
The integral kernel is equal to the field of a dipole of unit 
strength, 



Q(r,r'^) 



/i 2z 2 -{y-y') 



47T [ z 2 + ( r _ r /)2] 



./\2l5/2 ' 



(5) 



The integral Eq. (|4]) with kernel Eq. ([5]) is divergent at 
r — > y' and z — >• 0. In a numerical procedure, the diver- 
gence can be handled in three ways: (i) by keeping a finite 
z during the calculation^ (ii) by working in the Fourier 
space jii (iii) by converting the integral to a matrix form 
and using the flux conservation to determine the diagonal 
elements ] 12 i 16 i 17 Here we use the third method. Since, for 
H a = 0, the total flux through the z = plane is zero, 
the kernel should have the property J d 2 r Q(r, r', 0) = 0. 
This yields 



— B z (v) =H a +g(v)C(v) 
Mo 



L 



dV ff(r') - 5 (r) 



47T 



r'|3 



(6) 



FIG. 2: The simulated flux density map in a strip of width 
2w, with a semicircular indentation of radius ro = 0.2w, in 
applied field H a = 0.3J C , n = 19, and ramped with a rate 
HoH a — poJc/wd. Note that A is not equal to ro- 



where the scalar function C is an integral over the area 
outside the superconductor 



C(r) 



dr' 



utside 47r|r-r / | 3 " 
For a uniform strip of width 2w it yields 



Cstrip(^) — 



W 



(7) 



(8) 



In addition, the indentation gives a contribution from the 
semicircle, which is calculated numerically from Eq. ([7j)- 

In the following we use an equidistant square grid and 
ascribe the same area s to each grid point. The discrete 
version of the kernel then acquires the form 17 




(9) 



where = 1 / 4tt | — r j 1 3 f° r * 3 an d Qu =0. All ele- 
ments of the discrete kernel Eq. ([9|) are nondivergent and 
the flux conservation, J d 2 r B(v) = 0, is guaranteed. Re- 
lating the magnetic field and the current by the Faraday's 
law, and using the inverted Biot-Savart law one obtains 
the dynamic equation for the local magnetization: 



9(r,t) 
where 



d 2 r'Q-\r,v f ) \fg(r',t)-H a (t)} , (10) 



fg = V- (pVg)/dfi . 



For discrete formulation of the problem the inverse kernel 
Q~ l is just the inverse of the matrix Eq. ([9]), hence the 
matrix must be calculated and inverted only once. 



III. RESULTS AND DISCUSSION 

a. Magnetic field and current The simulations were 
performed by ramping the applied field at a constant 
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FIG. 3: Evolution of flux penetration depth a as the ap- 
plied field is ramped with a constant rate poH a — poJ c /wd. 
Stronger flux creep, i.e. smaller n, leads to deeper penetra- 
tion. The Bean limit is n — 101. 



rate noH a = poJ c /wd, starting at zero field and a flux- 
free strip. The flux penetrates from the edges forming 
well-defined flux fronts that move towards the strip center 
as the applied field increases. Shown in Fig. [2] is a typical 
result of the flux density distribution presented as seen 
in a magneto-optical image, i.e., the image brightness 
represents the magnitude of the perpendicular magnetic 
field. The sample edge is seen as a bright line, i.e., the 
flux density is highest at the edge. 

Far from the indentation the flux penetration front is 
straight, and leaves a fraction a/w of the strip in the 
flux-free Meissner state, seen here as a black region. The 
penetration of this straight front versus applied field is 
shown in Fig. [3] for different values of the creep exponent. 
For large n the simulations approach the Bean-model 
result, 4 aBean = w/ cosh(7ri^ a / J c ), while for smaller n, 
i.e., stronger flux creep, the penetration is deeper, all as 
expected for a strip with straight edges. 

Near the indentation the flux penetration largely fol- 
lows the circular shape. At both sides of the indentation 
there are dark regions of reduced flux density. As pen- 
etration gets deeper these will become narrow d-lines, 
where the current stream lines make sharp turns ^ In the 
Bean limit n —> oo, the d-lines of semicircular indenta- 
tions have parabolic shape. With finite n the parabolic 
shape is only approximated. However, the main effect of 
the indentation is that it pushes magnetic field deeper 
into the sample. In order to quantify this we define 
the excess penetration A as the difference between the 
deepest penetration and the penetration far away from 
the indentation. Fig. [4] shows how A evolves with in- 
creasing H a . Evidently, the excess penetration is not 
equal to the indentation radius, ro, as in the case of the 
bulk Bean model^^ Moreover, A turns out to be field- 
dependent. Initially, A increases, then reaches a maxi- 
mum followed by a decrease at larger H a . This surpris- 
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FIG. 4: Evolution of the indentation-induced excess pene- 
tration, A, as a function of applied field. The two panels 
correspond to different indentation radii, ro/w =0.1 and 0.2, 
respectively; poH a — poJ c /wd. 



ing non-monotonous behaviour is supported by magneto- 
optical measurements of the flux penetration in a uniform 
YBa2Cu30 :E film containing an edge defect, see Fig. [5l 
The film was shaped as a strip of half- width w = 0.4 mm, 
and the figure shows the flux distribution at 25 K for 3 
different applied fields. In (a) the field was very small, 
PoH a = 3 mT, creating negligible penetration so that 
the actual shape of the defect appears in the image as 
the bright "bay area" inside the strip. In this state the 
excess penetration is equal to the depth of the defect, and 
measures A = 80 pm. In (b) and (c) the applied field is 
17 mT and 36 mT, respectively, and the corresponding 
excess penetration is A = 115 pm and 100 pm. This 
gives for A/w = 0.20,0.29 and 0.25, demonstrating an 
excess penetration that exceeds the depth of the inden- 
tation by nearly 40 %, in very good agreement with the 
Bean model results plotted in Fig. [H 

The Fig. [4] includes the behaviour of A/w for two dif- 
ferent ro/w. Comparing the two panels we see that larger 
indentations produce a larger A. However, the relative 
excess penetration, A/ro is larger for the smaller indenta- 
tion. The excess penetration can exceed the indentation 
depth by almost 50 % for ro = O.lw and large values of 



FIG. 5: Magneto-optical images of flux penetration into an 
YBa2Cu3 0cc strip with a defect at the edge. Only the lower 
half of the strip is shown. In (a), (b) and (c) the applied fields 
was 3, 17 and 36 mT, respectively. The excess flux penetra- 
tion, A, is maximal at the intermediate field, in agreement 
with simulations. 



FIG. 6: Simulated flux distributions (top) and current stream 
lines (bottom) in an increasing applied field, where ro = O.lw 
and the other parameters the same as in Fig. [2] From left, the 
values of H a are 0.05J C , 0.2J C , and 0.4J C with corresponding 
values of a 0.98u>, O.Sw. and 0.5w. 



n. For smaller values of the creep exponent one always 
finds smaller A, implying that creep tends to smoothen 
perturbations in the flux front 4 

Our results demonstrate that an indentation in a thin 
film affects the flux distribution in a stronger and more 
complex way than it does in bulk superconductors. This 
must be due to the non-local electrodynamics of thin 
films, and in particular due to the presence of Meissner 
currents in the flux free regions. These Meissner cur- 
rents do not make the same sharp turns as the critical 
currents in the flux penetrated region, see Fig. [6] and also 
Refs-imi. As a result, the Meissner currents concentrate 
in front of the indentation where their density reaches j c 
and hence leads to even deeper flux penetration. This is 
why the flux front near the indentation advances faster 
than in the rest of the film. This accelerated advance- 
ment eventually terminates when the penetration depth 
becomes comparable to the strip halfwidth. The reason is 
simply that in the limit of full penetration all flux fronts 
reach the middle of the strip and hence A — > 0. 

b. Electric field The Lorentz force pushing magnetic 
flux is directed perpendicular to the local current density. 
Even a small indentation distorts the current stream lines 
over a large area, and hence significantly modifies the 
trajectories of flux motion. In particular, all the flux ar- 
riving to the fan-shaped region rooted at the indentation 
must have entered the sample through this indentation, 
see Fig. It creates a dramatic local enhancement of 
electric field since E is a direct measure of the intensity 
of flux traffic. 

Analytical solution for the electric field distribution 
around an indentation in thin films is not available. 
Therefore the results obtained for the case of a slab are 
often utilized as approximations also for films ^ 14 i 2Q We 
will now analyze to what extent such estimates are valid 
by comparing them with our simulation results for a strip. 

In the fan-shaped region that originates from the semi- 
circular indentation, the electric field can be found by 



solving the Maxwell equation V x E = — B in cylindrical 
coordinates. Since the evolution of ^-distribution is usu- 
ally not very far from the Bean model, one can assume 
B = jjLoH a , which leads to the solution^ 



(w - a + r ) 



w 



(w- \x\) 



(11) 



for \x\ > a — ro and zero for \x\ < a — ro. Far away 
from the indentation the solution of the same equation 
in cartesian coordinates, d x E = — fioH a , is 



E (x) = noH a (\x\ - a) 



(12) 



for \x\ > a and zero for \x\ < a. Note that the width 
w enters Eq. ([TTj) only because of the specific choice of 
the x-coordinate, where the edge is located at x = w. 
Replacement x — > x + w removes the dependence. 

Figure [71 compares Eq (x) and E\ (x) with the simu- 
lated electric field profiles. The quantitative agreement 
is poor, though the shape of profiles (both across the 
indentation and away from it) is fairly well reproduced, 
in agreement with Ref. 0. The expected enhancement 
of E due to indentation is also obvious. The formulas 
above predict the relative enhancement for the peak val- 
ues ^;( max ) ^( max ) — ^ w _ a )/2ro + 1 for a bulk sample. 
One can see from the plot that the effect of indentation 
is even stronger for thin films: the ratio ^( max ) /^( max ) [ s 
slightly higher and the excess penetration is larger (the 
flux front here corresponds to the point where E(x) = 0). 

A locally enhanced electric field near edge indentations 
and hence enhanced Joule heating is predicted to facili- 
tate nucleation of a thermal instabilit y 14 i 2Q The instabil- 
ity in thin superconductors is usually observed in form of 
macroscopic dendritic flux avalanches 2 ^ or macroscopic 
uniform flux jumps^ 2 -. However, a third scenario is also 
possible when a series of microscopic flux avalanches re- 
peatedly take place in the same region, each leading to 
a small advancement of the flux front ^ It creates an 
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FIG. 7: Electric field profiles across the indentation and far 
away from it. Solid lines: Eas. II II and 1121 for a bulk supercon- 
ductor. Dashed lines: simulations for a thin-film strip; the 
parameters are the same as for Fig. [2] except that ro = O.lw 
and H a /Jc = 0.25. E c = poJ c /d. A strong field enhancement 
near the indentation is clearly seen. 

additional front distortion since the avalanches are ex- 
pected to be larger and occur more frequently at the 
indentation, where the local E is maximal. Experimen- 
tally the individual avalanches can be very small, and 
hence it is not easy to determine whether the thermal 
effects contribute to an observed front distortion. To 
identify the penetration mechanism one can compare the 
observed flux profiles with the simulations. The maximal 
excess penetration due to non-thermal effects is found to 
be 150 % of the indentation radius for our parameters. 
Consequently, when the observed excess penetration is 
larger, the flux penetration probably occurs via thermal 
micro- avalanches . 



IV. CONCLUSIONS 

We have numerically solved the Maxwell equations to 
describe flux penetration into a thin superconducting 
strip with an edge-indentation and analyzed the time 
evolution of flux front in an increasing applied field, 
H a . The excess penetration, A, due to the indentation 
is not equal to the indentation radius, ro, in contrast 
to the well-known case of a bulk superconductor in the 
Bean model. Three different mechanisms that influence 



the excess penetration were analyzed, (i) The nonlocal 
electrodynamics in films leads to a characteristic A(H a ) 
dependence with a smooth peak. The ratio A/ro at the 
peak equals 1.5 when ro is 0.1 of the strip half- width and 
becomes even larger for smaller ro- (ii) The flux creep 
always tends to smoothen the flux front and decrease 
the excess penetration, (iii) Thermal flux avalanches 
are more likely to occur at the indentation, which can 
increase the apparent front distortion. Our results 
can be very helpful in order to identify which of these 
three mechanisms is the dominant one in a concrete 
experiment. 
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APPENDIX A: NUMERICAL DETAILS 

The simulations are carried out on an equidistant 
square grid with N x N points, x m = w(2m + i)/N — w 
and y n w(2n + 1)/N - w, for < m, n < N. The 
system has two symmetries that must be incorporated 
in the kernel: first, the periodic boundary, which means 
that we must add a mirror strips at x < — w and x > w. 
Second, the symmetry around x = 0. The latter means 
that we can work with half the kernel^ The simula- 
tions use a grid size of N = 100, which means that a 
5000 x 5000 matrix must be put in memory and inverted. 
The memory consumption is the main limiting factor of 
the simulations. The kernel is stable, so there is no need 
for additional smoothening. For most exponents a pure 
power law is used, but for the Bean limit, n = 101, a 
cutoff on the resistivity p < p ma x was necessary to en- 
sure stability. The flux front position was determined at 
every time step and then smoothened as a function of 
time. It allows the front position to be determined with 
an accuracy much better than the distance between two 
grid points. 
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